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ABSTRACT 


The  set  of  linearized  magnetohydrodynamic  equations  is 
reduced  to  the  reflectivity  equation  for  the  ccmpressional 
magnetic  perturbations  in  the  framework  of  the  Radoski 
model.  It  Is  shown  that  the  reflection  coeficient  is  a 
function  of  the  curvature  of  the  magnetic  field  lines,  the 
inhomogeneities  of  the  magnetic  field,  and  the 
inhomogeneities  of  the  Alfven  velocity.  An  interesting 
property  of  the  reflectivity  equation  is  that,  near  Alfven 
resonant  magnetic  force  lines,  this  equation  is  reduced  to 
the  curvature  free  Budden  equation.  We  show  that  near 
Alfven  resonances  the  curvature  does  not  play  a  significant 
role  and  3udden'3  assymtotics  In  time  can  be  applied  to  the 
wave  field  near  the  magnetic  force  lines  where  the  Alfven 
dispersion  relation  holds. 


1 


1.  PERSONNEL: 

The  following  specialists  are  involved  in  space  physics 
research  in  the  Geophysics  Program  of  the  Henry  Krumb  School 
of  Mines: 

1.  Philip  Carrion:  Wave  propagation  in  the  Earth's 

magnetosphere  and  Alfven  resonances. 

2.  Akira  Hasegawa:  Waves  and  instabilities  in  space 

plasma . 

3.  Waldo  Patton:  Numerical  approach  to  solving  the  MHD 

and  kinetic  equations.  Satellite  data 
processing . 

4.  Manju  Prakash:  Wave-particle  interactions. 

Professor  P.  Gross  who  is  an  expert  in  different  aspects  of 

plasma  theory  and  fluid  dynamics  serves  as  a  consultant  to 
the  project. 


2.  ACCOMPLISHMENTS: 

During  the  past  year  the  following  topics  were 
analyzed : 

1 .  Alfven  resonances  in  a  cold  plasma  magnetized  by  a 
curved  magnetic  field. 

Main  results: 

It  was  shown  that  a  set  of  linearized  MHD  equations  can 
be  reduced  to  the  reflectivity  equation  for  compressional 
magnetic  perturbation.  The  iefiection  coefficients  show 


that 


reflected  or  scattered  by  the 


the 


waves 


are 
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inhomogene  It ies  of  the  ambient  magnetic  field,  by  the 
curvature  of  the  magnetic  field,  and  by  the  inhomogeneities 
of  the  Alfven  velocity.  It  was  shown  that  near  the  Alfven 
resonant  field  line  where  the  Alfven  dispersion  relation 
holds,  the  reflectivity  equation  is  reduced  to  the  Budden 
curvature  free  equation  (Appendix  1)  and  that  Budden 
solutions  are  valid  near  the  Alfven  resonant  lines. 

2 .  Observation  and  analysis  of  Pc  3  -  Pc  4  pulsations  using 
the  geostationary  satellite  GOES-6. 

Main  results: 

It  was  shown  that  the  polarization  of  Pc  3  pulsations 
during  October  1984  (the  month  of  GOES-6  observations)  are 
polarized  linearly  in  the  azimuthal  direction.  Whereas  P°  4 
are  typically  elliptically  polarized,  our  study  of  wave 
characteristics  observed  at  the  spacecraft  indicates  that  Pc 
3  pulsations  can  be  caused  by  a  Kelvin-Helmholtz  instability 
occurring  at  the  magnetopause  whereas  most  observed  Pc  4 
pulsations  could  not  be  produced  by  this  type  of  instability 
and  indicates  that  an  internal  mechanism  generates  Pc  4 
pulsations  (Appendix  2). 

3.  Wave-particle  interaction  in  space  plasma 

Main  results: 

We  have  studied  the  response  of  charged  particle  flux 
in  ultra-low-  frequency  (UFL)  electromagnetic  field 
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superimposed  on  the  reference  dipole  field.  We  have 
investigated  both  the  linear  and  the  nonlinear  behavior  of 
particle  flux  by  tracing  numerically  the  motion  of  a  test 
particle  in  phase  space.  The  linear  response  has  been 
compared  with  that  obtained  analytically  from  the  kinetic 
equations.  The  calculated  particle  fluxes  were  compared 
with  the  observed  behavior  in  the  presence  of  geomagnetic 
pulsations.  Based  on  these  comparisons  we  inferred  wave 
characteristics  such  as  the  azimuthal  wave  number  "m" ,  the 
parallel  wave  vector  "k",  and  wave  localization  in  the 
radial  direction. 

Despite  many  experimental  efforts  to  infer  the 
properties  of  pulsations  from  particle  flux  oscillations, 
there  have  been  very  few  efforts  to  deduce  them 
theoreticaly.  Therefore  our  calculations  have  added  a  new 
dimension  to  the  theoretical  understanding  of  the 
phenomenon . 

In  addition,  our  calculations  were  used  to  study 
particle  diffusion.  We  have  also  studied  particle  diffusion 
arising  from  the  resonance  interaction  between  wave  and 
particle  motion  in  the  dipole  field.  Diffusion 
characteristics  are  considerably  different  from  those  of 
quasilinear  diffusion  because  of  chaotic  particle  behavior 
and  nonlinear  effects. 
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The  following  publications  have  been  sponsored  by  this  grant: 

1.  P.  Carrion,  "Alfven  waves  in  a  cold  plasma  with  curved 
magnetic  fields",  Phys.  Fluids  (in  press) 

2.  P.  Carrion,  "Compressional  and  toroidal  resonances  in 
the  magnetosphere",  (in  press) 

3.  P.  Carrion,  A.  Hasegawa,  L.  Lanzerotti,  "Study  of  Alfven 
standing  waves  at  a  geostationary  orbit",  (in  press) 
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3.  M.  Prakash,  A.  Hasegawa,  P.  Carrion,  ''Particle  Dynamics 
‘  in  Ultra  Low  Frequency  Field",  (Appendix  3) 
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ALFVEN  WAVES  IN  A  COLD  PLASMA 
WITH  CURVED  MAGNETIC  FIELDS 

Philip  M.  Carrion 

828  Selley  W.  Mudd  Building 
Henry  Krumb  School  of  Mines 
Columbia  University 
New  York,  NY  10027 


ABSTRACT 

The  set  of  linearized  magnetohydrodynamic  equations  is 
reduced  to  the  reflectivity  equation  for  the  compressional 
magnetic  perturbations.  It  is  shown  that  scattering  is 
caused  by  the  curvature  of  the  magnetic  field  lines,  by  the 
inhomogeneities  of  the  magnetic  field,  and  by  the 
inhomogeneities  of  the  Alfven  velocity.  An  interesting 
property  is  that  near  Alfven  resonant  lines  the  equation  is 
reduced  to  the  curvature  free  Budden  equation. 


INTRODUCTION 


Let  us  consider  the  set  of  linearized  MHD  equations 
that  govern  the  propagation  of  small  perturbations  in  a 
perfectly  conductive  magnetized  fluid: 


P£tt-  "  V  P  +  1  [(b-V)B  +  (B  . V)b]  (1  ) 

bo 

b'  •  -  7x  E  (2) 

t 

E  +[v  x  B]«  0  (3) 

b  -  -  B(V-C  )  +  (B-V  )  ?  -  (£*V)  B  (14) 


where  p  is  the  density,  £  is  the  3~D  displacement  vector,  B 
and  b  are  the  unperturbed  and  the  perturbed  magnetic  fields, 
respectively,  v  is  the  particle  velocity,  E  is  the  electric 
field  and  p  is  the  total  pressure: 


We  consider  an  anbounded  medium  and  assume  that  the 
Sommerf eld  radiation  condition  at  infinity  is  satisfied. 
The  initial  conditions  express  the  causality  of  the 
electromagnetic  field  and  the  particle  displacement  and 
velocity  that  the  field,  displacement  and  velocity  are  all 
zero  for  time  less  then  zero. 
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If  the  plasma  is  incompressible  (  .  o)  and 

homogeneous,  it  is  easy  to  obtain  from  (1)-(5)  that 

V2p  -  0  (6) 

Applying  the  divergence  theorem  to  (6)  yields: 

p  -  const 

Av  -  0  (7) 

where  A  is  the  local  scalar  Alfven  operator: 

32 

A  -tuoP  al2  ~  CB  -7  )2]  (8) 


If  the  plasma  is  not  homegeneous,  equation  (7)  is  no 
longer  correct.  Hasegawa  and  Uberoi1  considered  an 
inhomogeneous  cold  plasma  in  which  the  magnetic  field  lines 
of  the  unperturbed  magnetic  field  are  straight  but  the 
magnetic  field  is  not  homogeneous.  They  showed  that  in  this 
case  the  MHD  plasma  permits  both  singular  modes  and  a 
continuous  spectrum.  No  dispersion  relation  exists. 
Bartson2  showed  that  in  the  case  of  inhomogeneous  plasma,' 
Alfven  oscillations  are  damped  as  the  inverse  of  time. 

Uberoi  3  showed  that  the  governing  equation  for  Alfven 
waves  in  incompressible  fluids  is  similar  to  that 


governing 
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the  electrostatic  oscillations  of  a  cold  plasma.  Hasegawa 
and  Chen  4  and  Hasegawa5  considered  different  aspects  of 
Alfven-wave  propagation  in  nonuniform  plasmas.  They  treated 
the  problem  of  damped  and  undamped  surface  modes  for 
continuous  and  discontinuous  density  profiles.  Radoski  and 
McClay6  considered  the  propagation  of  Alfven  waves  in  cold 
plasmas  magnetized  by  a  curved  magnetic  field.  They 
discussed  a  possible  method  of  avoiding  spatially  dependent 
eigenf requencies  which  are  not  acceptable;  they  introduced 
an  electric  field  parallel  to  the  ambient  magnetic  field. 
Radoski7  solved  the  initial  value  problem  for  noncollective 
Alfven  resonances  and  showed  that  with  time  the  energy  in 
the  radial  mode  is  converted  into  azimuthal  motion. 
Recently  the  transient  i n i t i al- boundary  value  problem  for 
plasmas  with  curved  magnetic  fields  has  been  treated  by 

Q  ^ 

Carrion  and  Allan  et  al  .  They  used  tt.e  Radoski  model  and 
studied  numerically  the  conversion  of  Alfven  modes. 

The  knowledge  of  the  behavior  of  Alfven  waves  appears 
to  be  imprtant  not  only  for  the  laboratory  plasma  but  for 
space  plasma  as  well  helping  to  better  understand  low 
frequency  wave  propagation  (Lanzerotti  et  al10,  Chen  and 
Hasegawa1} Radoski  12  ) . 

In  this  paper  I  use  the  Radoski  model  which 
incorporates  a  cold  plasma  magnetized  by  a  magnetic  field 
whose  force  lines  are  circles.  First,  I  reduce  a  set  of 
linearized  MHD  equations  ( 1 ) - (5 )  to  the  reflectivity 
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equation.  The  reflectivity  equation  contains  a  term  which 
is  proportional  to  the  reflection  coefficient  of  Alfven 
waves  in  the  single  scattering  approximation.  This 
reflection  coefficient  is  a  function  of  the  curvature  of  the 
ambient  magnetic  field  and  of  the  i nhomogene i t i es  of  the 
density  and  the  Alfven  velocity.  The  reflectivity  equation 
is  transformed  into  the  Schroedinger  equation  whose 
eigenvalues  determine  the  coupling  between  two  motions  (in 
axial  and  radial  directions).  The  scattering  potential  of 
the  Schroedinger  equation  is  studied  and  the  solution  at  the 
singular  points  of  the  Schroedinger  operator  is  analyzed. 

REFLECTIVITY  EQUATION  AND  THE  REFLECTION  COEFFICIENT 

We  assume  that  the  unperturbed  magnetic  field  is 
circular.  In  the  cylindrical  coordinate  system  (R,cp  ,  Z)  we 
can  write  the  magnetic  field  as: 

B  "  <°»  0) 

(4)  into  (1)  and  assuming  a  dependency  of  the 
time  and  space  of  the  form: 

e  i(n<p  +  Az+  u  t) 

where  n  is  an  integer  and  A  and  u  are  real  numbers  we 
obtain : 


Substituting 
wave  field  on 
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(u>'y0P 


n 


M0V(  T  )  +  F 


where  F  is; 


(9) 


F  -  (B*V  )  [(  S-V  )  B  +  B  (  )  ]  -  (b-  V  )B  (10) 

Equation  (9)  shows  the  coupling  between  the  surface  mode 
and  the  Alfven  shear  wave  (toroidal  mode)  whose  dispersion 
relation  is: 

u2y  p  =  n2B2/R2 
0 


In  a  uniform  plasma  the  toroidal  mode  is  decoupled  from  the 
other  compressional  modes:  the  magnetosonic  and  the  ion 

acoustic  waves.  We  assume  that  the  plasma  is  cold  and 
consider  the  coupling  between  the  toroidal  mode  and  the 
poloidal  (compressional)  mode.  We  assume  that  there  are 
only  radial  inhomogeneities  in  the  density  and  in  the 
magnetic  field.  Since  the  plasma  is  cold  we  neglect  the 
pressure  P  in  (1).  Projecting  equation  (1)  onto  the  radial 
and  the  axial  directions,  I  obtain: 

0V1/B'  Ir  +  °I\, 

05j  -  l/B(iA  b^) 


(12) 


-6- 


where 

a-  |j  1»(« 

and  0  is  the  Alfven  algebraic  operator: 

0  =  (d2/v2  -  n2/R2 

A 


(n=  0,  1,2....  ) 


(13) 


(1H) 


where  v  is  the  Alfven  velocity.  For  a  °  0  and  R 

A 

(curvature-free  magnetic  field)  equations  (11)-(12)  become 
similar  to  those  discussed  by  Kivelson  and  Southwood13  .  We 


see  that  the 

shear  Alfven  mode 

is  coupled  with 

the 

compress i onal 

(poloidal)  mode. 

b  represents 

the 

compressional 

magnetic  perturbation. 

Elimination  £ 

R 

and  £ 

one  finds  an  equation  for  b^ 


b<p  ”  [ltl(9  /  +b(J>(  a/R  +  a'-  a(ln0  ) '+  0  _a2>=o^1^') 


where  all  derivatives  are  with 
is  known  in  wave  theory  as  the  " 
the  second  term  on  the  left  hand 
be  the  reflection  coefficient 


respect  to  R.  Equation  (15) 
reflectivity  equation"  since 
side  can  be  considered  to 
in  the  single  scattering 
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approximatlon ,  provided  that  the  radial  derivatives  of  9  and 
B  are  smooth.  Let  us  analyze  the  singularities  of  equation 
(15).  I  exclude  the  point  R  -  0  which  can  be  thought  of  as 
a  magnetic  shell  of  infinite  cuvature.  Then  the  obvious 
singularity  is  0(R  )  =  0.  The  Taylor  expansion  of  0 

yields : 

30 

0(R)  -  0 (R  )  +  _ 

3R 

R=R 

This  means  that 


(R  ~  Rfl)  +  0(|R  -  R0J  )  (16) 


[ ln(0) ] '  =  1/ (R-Rq) 


Setting 


X  -  R  -  R0 

I  obtain: 


V  '  -(1/x  +  (S  x-  A2-  a/x  +  ct/Ro+  a'  )  =0  (17) 

where  B  =  0 '  j 

R=R  o 

I  assume  first  that  near  the  Alfven  resonance  line,  a 

is  very  small  (  a  ■*  0).  Then  equation  (17)  can  be  reduced 
to  Budden’s  equation: 

V  "  (1/x)bi  +V3*-A2>=°  (18) 


which  is  the  same  as  for  the  curvature-free  magnetic  field. 
Now  I  can  draw  and  interesting  conclusiton  that  if  the  wave 
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is  propagationg  from  +» 

in  the  radial  direction,  it  will  be  reflected  or 

scattered  by  the  radial  inhomogeneities  of  the  magnetic 

field  and  by  the  ingomogeneities ,  the  Alfven  velocity 

profile  (Alfven  impedance)  as  well  as  by  the  curvature  of 

the  magnetic  field.  As  the  curvature  increases  the 

reflected  wave  has  larger  amplitude.  However,  near  the 

Alfven  resonance  line  the  governing  equation  is  reduced  to 

the  curvature-free  Budden  equation.  This  indicates  that 

near  the  Alfven  resonant  shell  the  curvature  does  not  play 

an  important  role  and  the  time  limit  solutions  for  b  and 

(P 

are : 


-*■  0  faster  than  ~ 

-  0  as  4 

The  Budden  solutions  applied  to  magnetospher ic  resonances 
were  discussed  by  Kivelson  and  Southwood13  .  They  introduced 
the  idea  of  cavity  resonances  in  the  magnetosphere.  An 
interesting  property  of  the  Budden  solutions  is  that  the 
magnetic  shell  which  corresponds  to  the  Alfven  resonant 
condition  behaves  like  a  reflecting  interface  with  an  energy 
transformation  occuring  in  the  magnetic  shell  (Alfven 

resonator).  The  problem  is  that  the  Alfven  resonant 
condition : 


v. 


n 

R 


n 


1  >  2  ,  .  .  . 


(19) 
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implies  that  a  number  of  magnetic  shells  will  resonate.  In 
general  (19)  represents  a  trancendental  algebraic  equation 


jy  *  \  00 

1  _  RtJ  (20) 

'  n 

The  first  equation  is  an  arbitrary  curve,  whereas  the 
second  is  a  straight  line  whose  slope  is  .  This  means 

that  the  intersections  of  the  line  and  the  curve  provides 
the  solutions  R  »  Ro  where  Ro  is  the  positions  of  one  of 
the  resonating  shells.  Suppose  now  that  the  source  is  not 
monohromatic  but  is  bandlimited.  Its  spectrum  then  can  be 
presented  as: 


a  (  to  )  - 


¥  0 
0 


u  C[  uL,  u1{  ] 

elsewhere 


where  and  Wh  are  the  low-cut  off  frequency  and  the 
high-cut  off  frequency  respectively.  Then  the  resonating 
shells  will  not  by  infinitesimal  but  of  a  finite  width. 
Carrion1-^  discussed  this  property  for  different  sources. 

Finally  equation  (15)  can  be  transformed  into  the 
Schroedinger  equation.  This  can  be  important  for  different 
aspects  of  Alfven  wave  theory  and  especially  for  the 
estimation  of  the  unknown  parameters  of  the  plasma  using  the 
observed  scattering  or  transmitted  Alfven  waves  • 

■  In  order  to  do  so  I  introduce  a  new  variable  u: 


10 


V  _9 

T  RB 


u  =  b„ 


Substituting  (21)  into  (15)  yields: 


(21  ) 


+  u  V 


-  x2 


\c  u 


(22) 


where  V  is  the  Schroedinger  scattering  potential: 

,2 


dR 


V 


RB 

0 


-  a/R  -  a'+a(ln0)’  -  0 


(23) 


RB 

T 


The  scattering  potential  is  singular  at  the  Alfven 
resonating  shells.  From  the  theory  of  Schroedinger  operator 
it  follows  that  strong  reflection  and  mode  conversion  usually 
occur  only  at  the  locations  of  the  singularities  of  the 
scattering  potential. 

The  eigenvalues  of  the  Schroedinger  operator  are  A  ; 

they  determine  the  coupling  of  the  radial  and  the  axial  particle 
motions. 

For  I  A  I  >>  I  v  I  (weak  scattering  or  large  coupling)  a 

solution  for  b  is 
<P 


V 


0/ (RB)  exp(-A  R) 


Which  decays  exponentially  with  R.  For  small  A  the 
coupling  is  small.  This  indicates  that  there  is  a  finite 
eigenvalue  A:(0<A<co)  which  provides  the  "best"  coupling. 
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CONCLUSION 

We  have  shown  that  the  set  of  linearized  MHD  equations 
can  be  reduced  to  a  single  reflectivity  equation  for  the 
compressional  magnetic  perturbation.  These  reflection 
coefficient  in  the  single  scattering  approximation  is  a 
function  of  the  curvature  of  the  ambient  magnetic  field,  and 
of  the  inhomogeneities  of  the  Alfven  velocity  as  well  as  of 

the  density  of  the  medium.  An  interesting  property  of  this 

equation  is  that,  near  the  Alfven  resonant  line,  it  reduces 
the  curvature-free  equation.  The  reflectivity  equation  can 
be  transformed  to  the  Schroedinger  equation  whose 
eigenvalues  determine  the  coupling  between  the  toroidal  and 
the  poloidal  modes.  The  Schroedinger  scattering  potential 
is  singular  at  the  location  of  the  Alfven  resonating 

magnetic  shells  which  means  that  at  these  positions 
reflections  and  mode  conversions  occur. 
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ABSTRACT 

Low  frequency  HM  waves  (in  the  range  of  Pc3  -  Pc4)  were 
measured  on  board  the  geostationary  satellite  GOES-6  during  the 
month  of  October  1984.  Fourier  analysis  of  the  three  components 
of  the  measured  magnetic  field  on  the  day  side  shows  that  the  Pc3 
pulsations  are  mostly  polarized  in  the  azimuthal  direction 
whereas  Pc4  pulsations  show  elliptic  polarization.  The 
cross-correlation  analysis  demonstrates  that  the  correlation 
between  the  transverse  component  of  the  wave  field  and  the 
component  of  the  magnetic  field  directed  along  the  ambient  field 
is  very  weak.  This  means  that  in  the  equatorial  region  where  the 
GOES-6  is  located,  no  compressional  modes  were  observed  in  the 
Pc3~4  ranges  on  the  day  side.  A  theoretical  model  which  predicts 
the  polarization  of  these  pulsations  is  discussed. 
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INTRODUCTION 

Alfven  waves  on  the  day  side  of  the  •  magnetosphere  can 
be  excited  by  a  variety  of  sources.  These  include  external 
sources  such  as  waves  excited  at  the  magnetopause  by  the 
Kelvin-Helmholtz  instability  and  internal  sources  such  as 
different  magnetospher ic  plasma  instabilities. 

The  main  characteristics  of  Alfven  standing  waves  have  been 
described  by  several  authors  (  see  for  example  Lanzerotti  and 
Southwood  1979,  Southwood  and  Hughes  1983,  Orr  1984).  The  data 
we  analyze  were  observed  by  the  GOES-6  which  carries  on  board  the 
Space  Environment  Monitor  (SEM)  instrument  package.  The  SEM 
includes  a  magnetometer  and  an  energetic  particle  detector. 
GOES-6  travels  in  a  geosynchronous  orbit  (6.67  Earth  radii)  at 
roughly  100  West  Longitude. 

We  analyzed  the  three  components  of  the  perturbed  magnetic 
field  on  the  day  side  in  the  Pc3-4  micropulsation  ranges 
(Pc3“10~45  sec.,  Pc4-45~150  sec.)  No  compressional  waves  in 
these  ranges  on  the  day  side  were  found.  However  the  azimuthal 
and  radial  components  of  the  Pc3-4  waves  found  were  strongly 
correlated;  they  represent  Alfven  shear  standing  modes.  The 
observed  Pc3  micropulsations  with  the  periods  near  30  sec.  were 
all  azimuthally  polarized.  However  the  observed  Pc4  pulsations 
were  elllptically  polarized  with  their  radial  component 
predominant.  The  rotation  of  polarization  as  a  function  of  local 
time  is  discussed. 
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THEORY 

Let  us  consider  a  cold  plasma  magnetized  by  an  ambient 
magnetic  field  directed  along  the  z-axis.  Assuming  that  the 
effective  dielectric  constant,  e,  as  well  as  the  magnetic  field 
vary  only  in  the  x  direction,  the  momentum  equation  combined  with 
the  frozen-in  condition  yields  (see  Hasegawa  and  Chen,  1976): 

c$—  r  3* 

cix.  ^  +  ^  Woc“  "f*0  8o  5(jc-*q\ 

a 

where  B0  is  the  magnetic  flux  density,  Wx  is  the 
x-component  of  the  3-D  displacement  vector,  o(u>)  is  a  Fourier 
transform  of  the  time  function  (source  wavelet)  and  x0  is  the 

location  of  the  source.  The  effective  dielectric  constant  can  be 
presented  by: 

£  (jc)  =  f „  ( u>z-  UJ*) 

(2) 

where  w0  is  the  Alfven  resonant  frequency: 

loo  —  Xa 

and 

i3  the  Alfven  velocity  and  p0 
the  Alfven  resonant  frequency 


is  the  mass  density.  Near 
(  “•  “  kz  VA  )  the 


(1) 
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effective  dielectric  constant  tends  to  zero. 

Let  us  consider  the  following  case.  Suppose  that  the  source 
is  a  surface  current  located  at  x  -  0  which  is  the  position  of 
the  magnetopause.  Then  inside  the  magnetosphere  near  the  Alfven 
resonant  frequency  equation  (1)  becomes: 


U/^-f 


ci.DC 


£ 


o 


(3) 


Assuming  that  |ej  <<  k*  BJ 

which  is  a  very  good  approximation  near  the  Alfven  resonant 
frequency  where  the  equation  (3)  can  be  reduced  to: 


f 


d? 

cLxd 


(4) 


where  rP  _ 

is  the  gradient  of  the  effective  dielectric  constant  near 
the  resonant  frequency.  A  general  solution  for  W  is: 

7*- 


W*=~  c  I 


cJLx. 


(5) 


where  Q  is  a  complex  function  which  does  not  depend  on  x 

The  next  step  is  to 


calculate  the  y-component  of  the 


displacement.  Since 
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d.  cc 


(  C  ->  o) 


(<> 


equation  (4)  can  be  reduced  to 


/  <£■_ 


C  Ux 


*} 


w  =-  o 
* 


(7) 


Equation  (7)  has  the  solution 


(8) 


where  ^  is  a  complex  function  which  does  not  depend  on  x 
Formula  (8)  shows  that  W  ^  has  a  simple  pole  at  the  resonant 
frequency,  which  implies  that  near  the  resonant  frequency 


_  /*>  cl 


/H/  / 


M/ 


ef 


6J 


(9) 


and  the  polarization  of  HM  waves  is  strictly  azimuthal  (the 
y-axis  represents  the  azimuthal  direction,  whereas  the  x-axis 
corresponds  to  the  radial  inward  direction). 

Under  the  assumption  of  a  linear  dependence  of  the 
dielectric  constant  similar  results  were  obtained  by  Hasegawa  et 
ai*  (1983)‘  The  integration  in  (5)  smooths  the  singularity  and 


this 

non- 

the 


means  that  even  if  the  effective  dielectric  constant  varies 
linearly  in  the  radial  direction  the  azimuthal  component  of 
hydromagnetic  pulsations  will  predominate  and: 

| Wy 1  >>  |WX|  near  e  -  0 


DISCUSSION 

Figure  (1)  shows  typical  hydromagnetic  pulsations 
observed  as  a  function  of  time  by  the  GOES-6  satellite.  The 
upper  curve  displays  the  compress ional  mode  (the  magnetic  field 
variations  directed  perpendicular  to  the  orbit),  the  middle 
curve,  the  radial  component ; and  the  lower  curve,  the  azimuthal 
component . 

Figures  (2)  shows  the  power  spectra  of  all  components  of  the 
magnetic  field  variation.  This  figure  illustrates  the  fact  that 
in  the  range  of  Pc3  pulsations  the  polarization  of  the  pulsations 
is  azimuthal.  (The  upper  curve  is  the  power  spectrum  of  the 
azimuthal  component  (letter  N).  Figure  3  is  the  hodogram  of  the 
Pc3  pulsation  presented  by  Figure  1  .  It  again  demonstrates  the 
azimuthal  polarization  of  the  Pc  3  pulsation.  During  the  entire 
month  of  October  1984  all  day  side  Pc  3  pulsations  were  polarized 
in  the  azimuthal  direction.  Our  observations  show  that  the 
observed  Pc3  pulsations  can  be  generated  by  a  narrow  band  surface 
current  located  at  the  magne-  topouse  due  to  the  occurence  of  the 
Kelvin-Helmholtz  instability.  Fgure  4  depicts  a  typical  Pc4 
micropulsation  as  a  function  of  local  time.  Figure  5  is  the 
power  spectrum  of  all  components  of  the  observed  geomagnetic 
pulsation.  The  upper  curve  is  the  power  spectrum  of  the  radial 
component  of  the  magnetic  perturbation.  This  indicates,  that  the 
observed  micropulsation  is  radially  polarized.  Our  obr  vations 
indicate  that  more  likely  the  Pc4  micropulsr.tions  can  be  also 
generated  by  the  Kelvin-Helmholtz  instability  and  the  resonant 
magnetic  tube  is  relatively  far  away  from  the  point  of 
observation.  If  the  point  of  observation  is  removed  from  the 


resonant  tube  the  polarization  of  HM  waves  can  be  elliptical  or 
even  predominantly  radial.  Figure  6  is  a  hodogram  of  a 
micropulsation  (PcH)  corresponding  to  the  morning  hours.  Figure 
7  is  an  observed  Pc  ^  pulsation  correspnding  to  the  afternoon 
local  hours.  Our  observations  indicate  that  statistically  a 
switch  of  polarization  at  about  local  noon  occurs.  This  switch 
of  polarization  is  intrinsic  for  the  Kelvin-Helmholtz  type 
instabilities.  Similar  results  were  reported  by  Lanzerotti  et 
al. (1 981  ) . 
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APPENDIX  3 


Particle  Dynamics  in  Ultra  Low  Frequency  Field 

M.  PRAKASH  (Columbia  University,  New  York,  NY  10027) 
AKIRA  HASEGAWA  (AT&T  Bell  Laboratories,  Murray  Hill, 

NJ  07974)* 

P.  CARRION  (Columbia  University,  New  York,  NY  10027) 

We  studied  the  behavior  of  particle  flux  in  ultra  Low 
frequency  (ULF)  electromagnetic  fields  superimposed 
on  the  reference  dipole  field.  Choosing  a  model 
hamiltonian  for  the  perturbing  field  we  have  traced 
the  particle  motion  numerically  in  two  dimensional 
phase  space.  These  calculations  are  performed  for  a 
wide  range  of  initial  energies  of  the  particle  typical 
of  those  in  radiation  belts.  The  particle  velocity 
is  studied  at  various  times  as  a  function  of  the 
amplitude  and  the  frequency  of  the  electromagnetic 
wave.  The  results  highlight  the  effects  due  to 
nonlinearity  and  wave-particle  resonance.  Our  calcu¬ 
lations  will  lead  us  to  study  particle  diffusion  at 
Landau  resonance. 
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